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In this work we study the contribution of the isoscalar tensor coupling to the realization of pseudospin 
symmetry in nuclei. Using realistic values for the tensor coupling strength, we show that this 
coupling reduces noticeably the pseudospin splittings, especially for single-particle levels near the 
Fermi surface. By using an energy decomposition of the pseudospin energy splittings, we show that 
the changes in these splittings come by mainly through the changes induced in the lower radial 
wave function for the low-lying pseudospin partners, and by changes in the expectation value of 
the pseudospin-orbit coupling term for surface partners. This allows us to confirm the conclusion 
already reached in previous studies, namely that the pseudospin symmetry in nuclei is of a dynamical 
nature. 
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I. INTRODUCTION 



Pseudospin symmetry is a concept that appeared in nuclear physics more than 30 years ago 0, Q to account for 
the observation, in heavy nuclei, of the quasidegeneracy of orbitals with quantum numbers (n, I, j = I + 1/2) and 
(n — 1, 1 + 2, j = 1 + 3/2) (for fixed n and I). Such pairs of single-particle states are known as pseudospin partners. This 
doublet structure is related to the pseudo-orbital angular momentum and pseudospin quantum numbers 1 = 1 + 1 and 
s = s = 1/2 0, 0> HI j respectively. The former, as noted by Ginocchio 0, is just the orbital angular momentum 
of the lower component of the Dirac spinor. Pseudospin partners are doublets with the same I. For example, for the 
partners [nsx/2, (n — 1)^3/2], I — 1, and for [np 3 / 2 , (n — l)/s/2] one has I = 2. 

The existence of degenerate or quasi-degenerate pseudospin partners is connected to a SU(2) symmetry of the 
Dirac equation with only scalar S and vector V potentials such that V = —S, regardless of the particular shapes of 
these potentials Q . It happens that in the Relativistic Mean-Field Theories of Nuclei (RMF) 0, ^J, 0] , the sum 
S = S + V is small at the nuclear energy scale, so that this symmetry, known as pseudospin symmetry, provides a 
natural explanation of the existence of quasi-degenerate pseudospin partners in nuclei referred to before |7l ll2Lll3Lll4) . 
However, since in RMF theories S acts as binding potential for the nucleons, it is not possible to have exact pseudospin 
symmetry in nuclei. Further works have shown that the particular shape of E, not just its smallness, can affect the 
pseudospin energy splittings and also explain the isospin dependence of these splittings [LH UJl llfl . Il7j . Moreover, 
using an energy decomposition coming from the Schrodinger-like equation for the lower component of the Dirac spinor, 
it was shown that the observed pseudospin splitting arises from a cancellation of the several ene rgy components, and 
not primarily from the pseudospin-orbit term, which is proportional to the derivative of E [TR llfij. Altogether, this 
led us to conclude that, in nuclei, pseudospin symmetry is realized in a dynamical way. A similar conclusion was 
reached by Marcos et al. 0, 0] . 

The tensor coupling has been used in studies of nuclear properties with effective Lagrangians including RMF 
theories by Furnstahl et al. in Ref. [23] , and in the relativistic Hartree approach model, studied by Mao in Ref. [2l| . 
Those works assessed its influence on nuclear observables, namely the spin-orbit splitting of single-particle levels in 
nuclei, the result being that the tensor coupling, a higher order term in a relativistic expansion, increases significantly 
the spin-orbit coupling. This suggests that the tensor coupling could have a significant contribution to pseudospin 
splittings in nuclei as well. This contribution is expected to be particularly relevant for the levels near the Fermi 
surface, since the tensor coupling depends on the derivative of a vector potential, which has a peak near the Fermi 
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surface for typical nuclear mean- field vector potentials. 

The tensor coupling has also been used as a natural way to introduce the harmonic oscillator in a relativistic (Dirac) 
formalism. In a recent paper, it was shown that the harmonic oscillator with scalar and vector potentials can exhibit 
an exact pseudospin symmetry |'2'2t |23|. When this symmetry is broken (S ^ 0), the breaking term is quite large, 
manifesting its nonperturbative behavior. However, if a tensor coupling is introduced, the form of harmonic-oscillator 
potential can still be maintained with £ = 0, but the pseudospin symmetry is broken perturbatively |24j . 

The tensor interaction has also been considered in order to explain how the spin-orbit term can be small for A- 
nucleus and large in the nucleon- nucleus case |2^. It is assumed that in the strange sector (case of A) the tensor 
coupling is large and the spin-orbit term obtained from this interaction can cancel in part the contribution coming 
from the scalar and vector interactions. This result shows that the tensor interaction can change strongly the spin- 
orbit term. In this spirit, we want to investigate if this interaction can also affect the pseudospin-orbit term (spin-orbit 
of the lower component) 0, @, IT^ . We address this problem by performing a mean-held calculation for the neutron 
levels of 208 Pb, using mean-field Lorentz vector and scalar potentials with a Woods-Saxon shape. These potentials 
were used in the previous works which revealed the dynamical nature of pseudospin symmetry. In this work we will 
perform a similar calculation, including now a tensor coupling term, and using again an energy decomposition similar 
to the one used in [l5j. We find that the tensor coupling potential has a noticeable effect on the pseudospin splittings. 
We study in particular detail the changes of the radial wave functions for the pseudospin partners and the differences 
between low lying and near to the Fermi surface pseudospin doublets. 

This paper is organized as follows. In Sec. Ill Al we present the Lagrangian for fermion fields coupled to external 
scalar, vector and tensor fields and obtain the corresponding single-particle Dirac Hamiltonian. The Dirac equations 
of motion are obtained in Sec. Ill Ul with emphasis on the second-order differential equations for the upper and 
lower components of the Dirac spinor. In Sec. Ill 01 we perform the energy decomposition based on the second-order 
differential equation for the lower component of the spinor, which will allow us to analyze the contribution of the 
radial tensor potential U for the pseudospin energy splittings. The results of the calculation using mean-fields with 
Woods-Saxon shape for S and V radial potentials (thereby also fixing the tensor potential) are presented in Sec. IIIII 
together with a discussion of the effects of the tensor potential on the neutron pseudospin partners, both for deep 
levels and levels near to the Fermi surface. The wavefunctions for the radial lower components of those levels are 
also plotted and their influence on the pseudospin splittings are discussed. Finally, our conclusions are summarized 
in Sec. |TV| 



II. DIRAC EQUATION WITH ISOSCALAR TENSOR COUPLING 



A. Dirac Hamiltonian 



Using the conventions of G. Mao [21|, the nucleon- meson Lagrangian density of a nuclear mean-field theory with 
nucleons interacting with cr, to and p mesons, in which a tensor (derivative) coupling is included, reads (h = c = 1) 



In this Lagrangian, u>^ v = d^u> v — d^u)^, p pv — d^p^ — d^p^ and u^ v = i/2[7 /1 , 7"]. Here the only vector mesons we 
are going to consider are the isoscalar to mesons. Furthermore, in a mean-field theory, the meson fields are static and 
only the time-like component is considered, i.e., we have 

cr = cr(r) lo^ = w(r)(5 M o ■ 

The Hamiltonian density is then given by 

H = 7-^ — rdo* - £ = *7°do* - £ = -iWa. V* + + g^a^ + g^V lu - i—^ (3a. Vw * . (2) 

d(do^) 2M 

The corresponding single-particle (Dirac) Hamiltonian is 

H = -i a. V + [3(M + g a a) + g u u - i ^Pa. Vw . 

ZM 

If we now define the scalar S and vector V potentials as, respectively, S — g a a and V = g U j^ 1 this last equation reads 



H = -ia.V + f3{M + S) + V - i [3 a U, 



(3) 
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where U — f w /(2M) Vw. Moreover, if the field uj(r) is just a function of the radial coordinate r, this Hamiltonian 
becomes 

H = -ia.V + /3(M + S)+V -i/3a-rU, (4) 

where U is the radial function 

In the remainder of the paper we will use the notation /„ = f^/g^ for the sake of simplicity. 

B. Equations of motion 

The Dirac equation for nucleons with tensor coupling is written as 

H^ = £^, (6) 

where H given by Eq. lj4"|l. 

It is instructive to decompose this equation into two second-order equations for the upper and lower components 
of the spinor iff, but retaining their spinor structure. To this end, we use the projectors P± = (I ± /3)/2 applied to 
i.e., define the spinors ^± = P± x &. Applying P± to the left of the Dirac equation © we get 

a ■ p + (M + S + V)S> + -ia-rUV- = £ (7) 
a + (-M -S + V)V- +ia-rUV + , (8) 

or still, defining E = V + S, A = V - S, and E = £ - M, 

{a-p-ia-rU)^- = (E — £)*+ (9) 
{a-p + ia-rU)^+ = (S + 2M-A)*_. (10) 

Using the formulas in the Appendix we get finally 

n 2 *_ + U7 2 - {/' - 2 h - 4(7 + 2 

F V r # - £/ £ - E dr \ E-EJ r 

= (E-T,)(E + 2M- A)#_ (11) 

2 T frrl rr, n U A' U \ A' 9*+ / r „ A' \ L ■ S T 

n 2 *+ + N7 2 + [/' + 2 — H f+ — — - + \ AU + 2 

y + \ r E + 2M -AJ + E + 2M - A dr \ E + 2M-AJ r + 

= (E - £)(£ + 2M - A)*+. (12) 

The terms with L ■ S for the upper and lower components are spin-orbit and pseudospin orbit coupling terms, 
respectively. We can note immediately that the tensor potential contributes to both. A full analysis will be done in 
the next section. 

If S, V and U are radial functions, then the general solution of Eq. (jfJJl is 



*«m(r) 



( \ 

r 

V — — y-«m(r) J 



(13) 



Here n is the quantum number related to the total angular momentum j and orbital momentum I by 

[ -0 + 1) =- (7 + 1/2), .7 = Z + 1/2 
K = S ■ ( 14 

1 Z = + (7 + 1/2), j = i - 1/2 



4 



The spinor spherical harmonics y Km result from the coupling of the two-dimensional spinors to the eigenstates of 
orbital angular momentum and form a complete orthonormal set. Through the relations 



J 



1 
2 

1 / K 

2VR 



(15) 
(16) 



one sees that, if the upper component of the spinor in Eq. I|13|) has a orbital quantum number I, the lower component 
(which has quantum number — k) must have a orbital angular momentum I = I — k/\k\. This quantum number has 
been associated with the pseudospin symmetry 0. 

Using the property er • ry Km = — 3^- K m, Eqs. © and ifTUI) reduce to a set of two coupled first-order ordinary 
differential equations for the radial upper and lower components g K and f K , namely, 



— + --U{r) 
dr r 



g K (r) = [£ + 2M-A(r)]/»(r), 



(17) 



dr 



U(r) 



Mr) 



-[E-E(r)]g K (r). 



(18) 



Similarly, from Eqs. (|11|) and l|12[) we arrive at the following second order differential equations for g K and f K 
k(k+1) A' 



dr 2 



E + 2M- A(r) 



d k TT , . 
— + - - U(r) 
dr r 



2 K^-U'{r)~U 2 {r)\g K (r) 



-\E -Y,(r)][E + 2M - A{r)]g K {r) 



(19) 



dr 2 



k(k — 1) 



£' 



E-S(r) 



C/(r) 



2 K ^ + {/'(r)-(7 2 (rU/ f; (r) 



= -[£-£(r)p + 2M-A(r)]/ K (r). 



(20) 



These two equations show explicitly the new terms that depend on U(r) and are originated by the tensor interaction. 
In particular, the term 2/t{7 (r)/r, which is the same for the upper and lower component, is the modification in the 
spin-orbit and pseudospin-orbit terms, respectively, generated by the tensor interaction. 

C. Energy decomposition and sum rule 

The terms in Eq. <|ll|l with denominator E — £, which have a singularity in E = S, fulfill a sum rule coming from 
Eq. l|A.2(i . If one divides each member of that equation by 2M* = E + 2M — A, left multiply them by and 
integrate, one gets 



T i a. ■ r . 

T £ d 3 r 



2M* 



~E-Z2M* dr 
P l L^_d 3 r-F /V — -^-^_d 3 r 



(21) 



E-E2M* J ~E-EM* r 

where P stands for the principal value of the integral. In terms of the radial functions g K and f K the sum rule reads 

p0 ° £' 1 - K 



£' f 00 E' 1 //_\' /"» £' U , 



■ 2M* 



E-E 2M* 



E - E2M 



E - E 2M< 



74 dr. (22) 



This sum rule can be used to check the numerical results. 

The energy decomposition of Eq. (jl 1|> can be performed by dividing it by 2M* and computing its expectation value 
for the spinor yielding 



( 2lF ) + {Vu) + {v * u) + ^ Da ™ in > + < Fpso > + ^ = E > 



(23) 
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where 



Vd a r wm 



(O) 



i jyu 

2M* E-Y, 

1 E' d 
2M* fi-E* 



V.SO = ^(-2U + -^-)^ (24) 



M* \ E-Y, 



\[rt ^3 r 



For the terms with E — E in denominator the integral is taken in the principal value sense. 



III. TENSOR COUPLING WITH MEAN-FIELD WOODS-SAXON POTENTIALS 

As stated above, the aim of this paper is to study the effect of the tensor coupling on the pseudospin splitting 
in nuclei. In previous works, in which we studied pseudospin symmetry in nuclei |14| . we solved numerically the 
Dirac equation with central mean-field potentials with Woods-Saxon shapes. Although these potentials are not full 
self-consistent relativistic potentials derived from meson fields, they are realistic enough to be applied to many nuclei. 
In this paper we follow the same approach, namely we consider the sum and difference potentials E and A to be of 
the general form 

p 

P{r) = l + cM(r-R)/a] ' (25) 
whereas the tensor potential U(r) is obtained by 

UM = Ay' = A S ' + A ' . (26) 
1 ' 2M 2M 2 K ' 

The depth, P , the radius (range), R, and the diffusivity, a, for E and A are fitted to reproduce the single-particle 
spectrum of 208 Pb [II Ej| . 

Using the general Woods-Saxon form in Eq. (|25|l for E and A, and with U given by Eq. I|26|l . we solved numerically 
the coupled first-order Dirac equations (|17|) and IjlH)). The single-particle energy levels are shown in Fig. 1, when /„ 
varies from to 1.3. This range of values is consistent to that one found in Ref. |2Qj for fittings of RMF and point 
coupling models to nucleon observables. 

One sees clearly from Fig. 1 that turning on the tensor coupling decreases the pseudospin splittings of the levels 
near the Fermi surface. For f v — 1.3 the pair [2/5/2 — 3^3/2] becomes almost degenerate while the [lin/2 — 2^9/2] 
doublet even reverses its order. Interestingly, the experimental energy values for these two pairs show this order 
reversal, which is not reproduced by the model calculations without tensor coupling. In Fig. ^ when we increase the 
tensor coupling f v , the energy for pseudospin partners with n < become deeper and those with n > become more 
unbounded. This systematics indicates that pseudospin symmetry is improved by the tensor interaction. 

In order to understand better why and how tensor coupling affects pseudospin splittings, we computed the contri- 
butions from the terms of the energy decomposition in Eq. i|23[) to the energy splittings for both low-lying and close 
to the Fermi surface pseudospin partners. 

In Table [I] we can see what is the contribution of all the terms of that energy decomposition for the two lowest 
neutron pseudospin partners ([2si/ 2 — W3/2] and [2p 3 / 2 — I/5/2]) an d f° r the t w0 topmost neutron pseudospin partners 
([3p3/2 — 2/5/2] and [2gg/2 — 1*11/2]) f° r tw° values of the tensor coupling strength: /„ = 0.0 (no tensor coupling) 
and f v — 1.3. As expected, since tensor interaction is an higher order interaction in the Lagrangian (derivative term) 
scaled by 1/M, the changes in the energy produced by the potential terms Vjj and Vs'U are small in comparison with 
the kinetic and potential terms, as we can see from Tab.|I] This table also shows that these terms, together with the 
pseudospin-orbit term, are significantly bigger for the surface levels than for the lower levels. This agrees with our 
expectations, referred before, that the effect of the tensor coupling is larger for the surface levels, since the potential 
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FIG. 1: Calculated neutron single-particle energy levels of the pseudospin partners [2/7/2 — I/19/2], [2/5/2 — 3p3/2] and [H11/2 — 
259/2] m 208 Pb. The left-most values are the corresponding values of the model Gl of Furnsthal et al. 26] and the experimental 
values ,2Tl ar e in the right-most column. In the middle column (WS) is our calculation with Woods-Saxon potentials for three 
values of the tensor coupling strength f v , corresponding to lines of different thickness. Pseudospin partners levels with k < 
are represented by full lines whereas those with k > have dashed lines. The Woods-Saxon parameters used to fit 208 Pb 
neutron energy levels are R = 7 fm, Aq = 650 MeV, £0 = —66 MeV and a = 0.6 fm. 



U is proportional to the derivative of the vector potential. On the other hand, changes induced by the tensor coupling 
in (Vs'u) and (Vbarwin) for surface levels are smaller. 

The smallness of the terms containing the tensor potential in regard to the kinetic and E potential terms is misleading 
concerning the effect in the pseudospin symmetry, since the changes of these last terms with f v are small, whereas the 
corresponding changes in (Vpso) and (Vu) can be significant for the surface levels, especially when compared with 
the energies of these levels. Note that the values of (Vpso) can have quite different values for different levels. These 
is due to the respective values of the k quantum number, as explained below. 

To have a better understanding of the influence of all these terms in the pseudospin splittings, we plot in Figs.[3Ja) 
and [21 b) the splittings and the differences of the terms in (|23[) for the pseudospin doublets [2si/2 — 1^3/2] and 
[2.99/2 — 1*11/2] as a function of /„. From these figures is clear that the decrease of the pseudospin energy splitting AE 
is much more pronounced for the surface doublet than for the low-lying one, and would be even more if we considered 
the relative energy variations. 

A more detailed analysis of the several contributions to the pseudospin energy splittings reveals that for the deep 
pseudospin doublet the contribution of the pseudospin-orbit potential Vpso almost does not change with /„, whereas 
the contribution of terms like (p 2 /2M*) and (E), which do not depend explicitly on the tensor potential U, is greater. 
This means that the main contribution to the change of pseudospin splitting with the strength of the tensor potential 
comes mainly via the change of the wave function induced by U. Furthermore, the energy splitting results for the 
most part from a cancellation of the A(E) and A(Vbarwin) contributions and also from A(p 2 /2M*) , the contribution 
from Vpso having a lesser role. This agrees with previous findings of similar studies of pseudospin splittings [l5l,llt|. 

As far as the upper doublet [2<?g/ 2 — H11/2] is concerned, much of the previous analysis still holds, except for the 
fact that the pseudospin-orbit potential is much stronger and also changes sensibly as /„ changes, being responsible 
for most of the pseudospin splitting. The reason for this, and in particular the fact that for values of f v greater than 
0.95 the splitting becomes negative, i.e., there is a level inversion, must be found in a more detailed analysis of the 
contribution of the tensor potential U for the Vpso potential. From Eq. I|24l) we see that, while U gives a positive 
contribution to this potential (since V < 0, so that, from Eq. Q, —U > 0), the effective contribution depends on 
the sign of E'/ (E — E) x L ■ S/r. Since it changes sign when r > r s , in which r s is the radius for which E = E, the 
net change of (Vpso) depends on how the wave function f K (r) behaves near r s for both pseudospin partners. What 
we found is that for these surface levels, the U contribution for the pseudospin-orbit potential 



Vpso- M .U r 



(27) 
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partners 


fv 


< P 2 /2Ar> 


(Vu) 


<vw> 


(^Darwin) 


<Vpso) 


(£> 


E 


2S1/2 


0.0 
1.3 


24.4396 
23.9037 


0.0000 
0.2351 


0.0000 

-0.0943 


-3.9527 
-3.6870 


-0.5852 
-0.5632 


-61.4644 
-61.9114 


-41.5627 
-42.1170 


14/2 


0.0 
1.3 


21.1032 
20.7075 


0.0000 
0.1459 


0.0000 
-0.0712 


-0.8106 
-0.4170 


0.0966 
0.1559 


-64.4159 
-64.6181 


-44.0266 
-44.1170 


2P3/2 


0.0 
1.3 


33.2950 
32.7195 


0.0000 
0.2465 


0.0000 
0.1223 


-2.7538 
-2.7319 


-1.6340 
-1.5835 


-60.2687 
-60.8762 


-31.3615 
-32.1033 


1/5/2 


0.0 
1.3 


28.5114 
28.2218 


0.0000 
0.2598 


0.0000 
-0.0021 


0.7303 
0.9522 


0.5384 
0.6780 


-64.5165 
-64.6667 


-34.7365 
-34.5571 


3P3/2 


0.0 
1.3 


52.9135 
52.6913 


0.0000 
0.5000 


0.0000 
0.1669 


0.8482 
0.9638 


-1.6013 
-1.9231 


-59.8025 
-60.2278 


-7.6420 
-7.8280 


2/5/2 


0.0 
1.3 


50.1727 
50.0072 


0.0000 
0.6256 


0.0000 
0.1427 


2.3696 
2.3831 


1.0254 
1.3089 


-62.2874 
-62.3270 


-8.7197 
-7.8595 


2ff9/2 


0.0 
1.3 


55.7666 
56.0706 


0.0000 
0.5622 


0.0000 
0.2855 


2.8295 
2.6487 


-6.4555 
-7.2110 


-53.3816 
-54.4705 


-1.2410 
-2.1140 


1*11/2 


0.0 
1.3 


51.2033 
51.6221 


0.0000 
0.8997 


0.0000 
0.1529 


3.1788 
2.9268 


3.3530 
4.1391 


-61.0308 
-61.2232 


-3.2958 
-1.4826 



TABLE I: Values of energies and terms in the decomposition 12311 for the pseudospin partners [2si/2 — 14/2], [2p3/2 — I/5/2], 
[3^3/2 — 2/5/2] and [2gg/ 2 — lin/2] for two values of f v . The energies and expectation values are given in MeV. 

is the dominant one. The reason why the contribution of this potential to AVpso for the pair [2g 9 / 2 — 1*11/2] is 
so big lies in the fact that it is proportional to (2L ■ S) = — (1 — k), which is equal to —6 and 5 for 2g 9 / 2 and 
1*11/2 respectively. Since U is negative, (Vp^ )(2gtg/ 2 ) — (^psoH-^n/z) is negative, and more so as f v increases, (see 
Fig. |2Ib)), thus explaining why pseudospin splitting decreases with the increase of /„. 

Note that, from Tabled one is able to see directly the different effects that the tensor coupling has on A(Vpso) 
in a low- lying ([2p 3 / 2 — I/5/2]) and a surface ([3^3/2 — 2/ 5 / 2 ]) pseudospin partner that have the same k values. The 
values of (Vpso) and their differences change significantly with f v for the surface partner. 

In Figs. Ota) andEfb) are represented E, Vu, the radial part of the Vpg Q potential [Eq. lj2"?j) ]. given by —U / (M* r), 
which is always positive as discussed before, and the lower radial wavefunctions f K for the pseudospin partners 
[2s!/ 2 — I4/2] and [2g 9 / 2 — 1*11/2] respectively, when f v = 1.2. Also plotted is the energy of the lowest lying level of 
each doublet for f v — 0, allowing to have a rough estimate of the value r s mentioned above by the intersection of its 
horizontal line with the £ potential curve. From these figures one sees clearly why the tensor potential U has a much 
larger effect for a surface pseudospin doublet. Indeed, the lower radial wavefunctions for these levels have a significant 
strength near the nucleus surface, such that (Vpg ) can have a significant value and have a sizeable contribution to 
the pseudospin splitting. 

IV. CONCLUSIONS 

In this paper we have assessed the importance of the isoscalar tensor coupling to pseudospin symmetry as is realized 
in heavy nuclei, using as an example the neutron pseudospin partners of 208 Pb, calculated within a relativistic theory 
with scalar and vector mean-fields parametrized with Woods-Saxon potential forms. By looking into the second-order 
equation for the lower component of the Dirac spinor, we obtained an expression for the pseudospin-orbit potential, 
showing its explicit dependence on the radial tensor potential. Since this potential comes from a derivative coupling 
(higher order term in the Lagrangian) the contribution for the energy from the potentials originated by that interaction 
are all scaled by 1/M (M is the nucleon mass) and, because of that, is very small in comparison with kinetic and 
potential terms. However, since the spin-orbit interaction (and pseudospin) is a term of the same order, the effect of 
the tensor coupling can be significant in this case. In fact, we conclude that the contribution from U to the pseudospin 
potential Vpso is the dominant one for the surface levels. 

We have shown that the surface pseudospin partners were the most affected by the tensor coupling, as was expected, 
and that this coupling reduces pseudospin splitting. This reduction can be significative to the point of inverting the 
level order, with the states with aligned spin (j = 1 + 1/2) having higher energy than the states with anti-aligned spin. 
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FIG. 2: The contributions from the terms of the energy decomposition (ll'MI for pseudospin doublets (a) [2s!/ 2 — ld3/2] and (b) 
[2ff9/2 — I211/2] when /„ varies from to 1.3. 



By analyzing the several contributions to the pseudospin energy splittings, we were able to confirm a conclusion 
reached in previous works, namely that the pseudospin symmetry is realized dynamically in nuclei, resulting from a 
cancellation of the various contributions, rather than just the one from the pseudospin-orbit potential. We also found 
a systematic change in the energy for the pseudospin partners when we increase the tensor coupling f v : states with 
aligned spin (k < 0) become deeper in opposite to anti-aligned ones that become more unbounded. This systematics 
allows us to conclude that pseudospin symmetry is improved by the tensor interaction. Finally, we were able to show 
for a surface pseudospin doublet the pseudospin-orbit potential, especially through its tensor potential part, gives a 
significant contribution for the change of the energy splitting. This finding is compatible with the dynamical character 
of the pseudospin symmetry. 
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FIG. 3: The potentials E, Vc/, —U/ (M* r), and the lower radial wave functions f K for the pseudospin partners (a) [2si/2 — ld3/2] 
and (b) [2<7g/ 2 — lin/2] when /„ = 1.2. The wave functions were normalized such that they could be plotted side by side with 
the potentials. 
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APPENDIX 



In this Appendix we present the derivation of equations (|llfl and (|12|) . 

If we define the operator O = ap + iarU, the second-order differential equations are obtained by applying OO* 
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and O^O to and *S? + , respectively. For we have 

00t^r_ = 0(£J-E)*+ = [o!-p(£J-E)]* + + (£;-S)C>* + 
= i a • f E' + (£ - E)(£ + 2M - A)*_ 

p 2 + C/ 2 -[/'-2^-4C/^-^^_ , (A.l) 

where it was used the fact that E, A and U are radial potentials, and primes denote derivatives with respect to r. S 
stands for the 4x4 spin matrix. On the other hand, from the equation O'ty- = (E — we obtain 

O f y>- = (oc-p-ia-rU)^- 

= a ■ f(f ■ p + if x p ■ E — iU)^f- 

■ TTT , L S T 

= tar — r-t/*_-2 

\ or r 

= (£?-£)*+ , 

which allows us to write 

0t\i) e' / L S 

ia-fE'* + =ia-fE'^^ = -^f^ + f/*_-2^*_ ) . (A.2) 



In the same way, for ^ + we have 



O t O* + = ia-rA'*_ + (£J-E)(£J + 2M-A)* + 

p 2 + (7 2 + C/' + 2-+4C/^)* + , (A.3) 



r r 



and, using equation 0$ + — (E + 2M — A)*-, we get 

* — A '*- = E + 2M-a (^-^- 2 ^)- 
Using the equations (|A.1|I and (|A.2|I for and (|A.3|I and (|A.4|I for W+, we get finally the equations lfTT|) and 
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